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1.  INTRODUCTION 


Passive  satellite  sounders  do  not  directly  measure  wind  velocity;  instead  they  derive  it 
indirectly  from  the  satellite  temperature  profile.  The  temperature  profile  is  computed 
from  radiances  (infrared  or  microwave)  emitted  by  volumes  of  atmosphere,  each  defined 
by  a  horizontal  resolution  area,  and  a  vertical  thickness  according  to  a  weighting  function 
centered  at  some  pressure  level.  The  temperature  profile  is  used  to  compute  geopotential 
heights  for  the  several  pressure  levels.  The  height  field  from  many  soundings  for  each 
pressure  level  is  used  to  compute  the  geostrophic  wind.  The  geostrophic  wind  is  an 
approximation  to  the  real  wind  that  assumes  frictionless  flow  parallel  to  "straight"  height 
contours  (along  a  great  circle),  where  the  pressure  force  balances  the  coriolis  force.  An 
alternative  method  involves  the  shear  of  the  geostrophic  wind,  commonly  referred  to  as 
the  "thermal  wind."  For  a  given  layer,  thermal  wind  may  be  derived  by  using  the  gra¬ 
dient  of  mean  layer  temperature  or  the  gradient  of  layer  thickness.  Since  the  thermal 
wind  is  really  a  shear,  the  wind  velocity  at  some  baseline  level  must  be  measured  or 
computed  (for  example,  geostrophic  or  measured  wind  at  800  hPa).  Figure  1  illustrates 
the  meaning  of  the  thermal  wind  through  a  simple  example.  In  the  example  the  assump¬ 
tion  is  that  the  geostrophic  wind  at  700  hPa  (1  hPa  =  1  mbar)  has  a  speed  and  direction 
of  20  ms'1  and  240°,  and  at  500  hPa  the  values  are  20  ms'1  and  300°.  The  thermal  wind 
speed  and  wind  direction  are  20  ms 1  and  360°,  respectively.  In  this  case  the  isotherms 
run  from  north  to  south,  with  the  warm  air  to  the  west.  Geostrophic  wind  and  thermal 
wind  are  defined  and  derived  in  standard  texts  (Haltiner  and  Martin  1957,  Holton  1979). 


Figure  1.  Illustration  of  the  thermal  wind  for  a  layer  from  700  to  500  hPa  (see  text). 
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This  report  presents  some  information  and  sample  calculations  on  the  types  of  errors  that 
may  be  expected  when  geostrophic  wind  or  thermal  wind  is  used  to  estimate  the  actual 
wind  from  data  gathered  by  passive  satellite  sensors.  The  assumption  is  that  no 
additional  data  are  available  (that  is,  only  satellite  sounding  data  from  sensors  of  the  type 
found  on  present  day  environmental  satellites).  This  report  does  not  treat  the 
improvements  that  should  occur  when  data  from  other  sources  of  data  (for  example, 
radar  profilers,  ground-based  radiometers,  unmanned  aerial  vehicles)  are  combined  with 
the  satellite  data. 


2.  GEOSTROPHIC  WIND 
2.1  Basic  Derivations 


The  equation  for  geostrophic  wind  along  a  constant  pressure  surface  is  given  by 


1  3$ 
fdn 


(1) 


where  f  is  20  sin  0  (0  sin  0  is  the  angular  velocity  of  the  earth  about  the  local  zenith 
at  latitude  0 ,  and  O  is  the  angular  velocity  at  the  poles),  n  is  the  horizontal  distance, 
<f>  is  geopotential  height,  and  V  is  the  geostrophic  wind  component  perpendicular  to  the 

n  direction.  In  general,  geopotential  height  is  4>  =  where  z  is  height  and  g  = 

gravitational  acceleration.  For  z  in  meters  $  *  0.98z. 


The  geostrophic  wind  also  may  be  computed  for  a  constant  level  (height)  surface  where 
1  dp  l  KT  - 

V  = - —  where  —  =  — -  (T  =  mean  virtual  temperature  along  the  distance  n  on 

*  p/a»  p  p 

the  surface,  p  =  mean  pressure  along  that  same  path,  and  R  is  the  gas  constant  for  dry 
air).  The  virtual  temperature  accounts  for  moisture  in  the  air,  which  having  a  lower  den¬ 
sity  than  dry  air  causes  Tv  >  T.  However,  at  temperatures  normally  found  above  4  or 
5  km  T  *  Tr  and  for  the  purposes  of  this  report  one  may  consider  T  *  Tv  at  satellite 

heights  >2  km  (generally  true  over  land  except  in  a  moist  summer  or  tropical  atmo¬ 
sphere).  This  second  form  of  the  geostrophic  wind  equation  can  be  used  to  check  for 

changes  caused  by  errors  in  p  or  Ty  (or  T). 


A  similar  type 
gradient  wind. 

equation  is  Vv 


of  "equilibrium"  wind  for  curved  parallel  contours  is  referred  to  as 
In  this  case,  the  tangential  acceleration  =  0.  The  scaler  version  of 


2 


-1  +  1  + 


where  r  =  R  tan  «  (R  =  radius  of  earth, 


the 

the 


a 
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=  angular  radius  of  the  equivalent  small  circle  for  r  as  seen  from  the  center  of  the  earth, 
and  r  =  radius  of  curvature.)  The  equation  may  be  rewritten  for  V  and  V  in  knots  and 
nautical  miles,  or  kilometers  per  hour  and  kilometers. 


V  -  180Grf  I  -1  +1+4 


r,  li 


3600 rft 


The  parameter  rf  may  be  computed  or  extracted  from  standard  tables  to  a  sufficient 
accuracy  (List  1984). 

Heights  of  pressure  levels  are  normally  computed  from  the  hydrostatic  equation 


where  the  subscript  o  refers  to  the  surface  or  bottom  of  a  layer.  As  before,  T  may  be 
substituted  for  Tr.  Here  the  mean  value  is  along  a  vertical  path  for  a  layer  of  atmo¬ 
sphere.  Computation  begins  at  the  surface  or  a  reference  p  level.  The  height  differences 
or  layer  thicknesses  (i  -  Zq  *  at)  are  "stacked"  upon  each  other  from  bottom  to  top. 


The  real  wind  may  be  defined  as  the  geostrophic  wind  plus  an  ageostrophic  component. 
The  horizontal  part  of  this  component  consists  of  a  pressure  tendency  term  and  a  term 
involving  the  change  in  spacing  of  the  contours  (or  isobars).  To  a  sufficiently  reasonable 

approximation  (Haltner  and  Martin  1957),  the  pressure  tendency  (—1  term  is  given  by 


where  s  consists  of  normal  and  tangential  distances.  The  term  for  the  change  in  contour 

ydV 

spacing  may  be  found  by  (V  -  Vg^  «  — where  V  is  the  real  wind.  To  a  fairly 

reasonable  approximation,  V  on  the  right  side  may  be  approximated  by  Vt,  the  mean 
value  along  the  distance  s  (in  this  case  along  a  path  equidistant  from  neighboring 
contours). 
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•  5(5)  (5) 

One  may  compute  values  that  roughly  estimate  "typical"  magnitudes  of  the  possible 
errors  in  wind  caused  by  the  geostrophic  approximation  and  errors  in  the  geostrophic 
wind  caused  by  errors  in  input  values  (for  example,  mean  layer  temperature).  As  a  first 
step,  wind  speed  only  is  considered. 


2.2  Errors  In  Vf  Caused  By  Mean  Profile  Temperature  Error 

The  geostrophic  wind  for  contours  at  100  geopotential  meters  (gpm)  in  units  of  knots 

001712 

(List  1984)  may  be  stated  as  V0  «  — - where  An  is  in  degrees  of  latitude.  At  45* 

*  A/t 

latitude  and  a  contour  interval  of  30  gpm,  and  An  expressed  in  kilometers, 


„  5.701  x  10»  .  . _ _  2.9387  x  10*  . 

V  m  -  m  kn  or  V  »  - in  ms 

♦  AM  • 


An 


a n 


(6) 


For  other  contour  intervals  multiply  by  the  ratio  of  "new  interval"/30  gpm.  For  a  typical 
midlatitude  atmosphere  (Jursa  1985),  for  the  surface  (assume  at  1000  hPa)  to  500  hPa 
layer,  the  possible  error  can  be  computed  from  net  temperature  profile  error  by  using 

equations  (3)  and  (6).  The  mean  layer  temperature  for  the  April  atmosphere  at  45CN 
(Jursa  1985)  is  used  to  compute  the  height  of  the  500  hPa  level. 


AZ 


-  287  K) 

9.81  ms'2 


(-  0.693) 


5373.6  m  *  5266.1  gpm 


Values  are  computed  for  differences  of  T  of  1  and  2  K,  or  T  =  264  and  263  K. 
a *  5353.3  m  *  52462  gpm 
aZjo  *  5333.1  m  *  5226.4  gpm 

For  the  case  of  a  difference  over  the  1000  to  700  hPa  layer  for  the  same  atmosphere, 

=  :2X7_r'm?s-2  *)(_  0.356675)  =  2797.0  m  *  2741.1  gpm 

9.81  ms'2 


a Z20  =  2786.6  m  *  2730.9  gpm 
a Zjfg  *  2776.1  m  *  2720.6  gpm 
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Similar  computations  may  be  made  for  other  layers  or  sublayers,  for  example,  in  table  1. 
We  now  assume  a  "true"  a$  of  30  gpm  between  soundings  200  km  apart.  The 
differences  in  T  errors  between  the  soundings  are  assumed  to  be  1  and  2  K  (for 

example,  the  first  s  jnding  has  a  T  error  of  2  K,  and  the  second  has  an  error  of  3  K, 
for  a  difference  in  error  of  1  K). 

TABLE  1.  DIFFERENCES  IN  GPM  COMPUTED  FOR  THE  LISTED 

LAYERS  FOR  THE  LISTED  DIFFERENCES  IN  T  ERROR 
BETWEEN  TWO  SOUNDINGS  200  KM  APART. 


T  Difference 

Layer  (hPa) 

(K) 

1000-900 

1000-800 

1000-700 

1000-800 

1000-500  | 

0 

0 

0 

0 

0 

0  1 

1 

3.0 

6.4 

10.2 

14.7 

19.9  | 

2 

6.1 

12.8 

20.5 

29.3 

39.7  | 

Equation  (6)  can  be  used  to  compute  the  Vg  values  for  the  "measured"  difference  (30 

gpm)  over  200  km,  and  Vg  arising  from  the  two  T  error  differences  (1  and  2  K),  for 
both  the  entire  1000  to  500  hPa  layer  and  the  several  sublayers.  Table  1  shows  the  gpm 
differences  and  table  2  shows  values  of  Vg . 

TABLE  2.  VALUE  OF  Vg  IN  ms1  COMPUTED  FOR  THE  LISTED  LAYERS 

AND  DIFFERENCES  IN  T  ERROR  BETWEEN  TWO  SOUNDINGS 
200  KM  APART,  a*  DIFFERENCES  WERE  ADDED  TO  THE 
"MEASURED"  30  GPM  TO  OBTAIN  THE  VALUES  SHOWN. 


T  Difference 

Layer  (hPa) 

(K) 

1000-900 

1000-800 

1000-700 

1000-800 

- - - -y 

1000-500 

0 

14.7 

14.7 

14.7 

14.7 

14.7 

1 

16.2 

17.8 

19.7 

21.9 

24.4 

2 

17.7 

21.0 

24.7 

29.0 

34.1 
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Table  2  shows  that  the  V  differences  are  about  9.7,  5.0,  and  1.5  ms~l  respectively  for 

the  layer  and,  for  example,  the  1000  to  700  hPa  and  1000  to  900  hPa  sublayers  for  a  T 
error  difference  of  1  K.  For  a  2-K  difference,  the  Vg  differences  are  about  19.4,  10.0, 

and  3.0  ms~l,  respectively. 

These  differences  in  V  are  not  small  and  may  even  result  in  a  180-degree  change  in 
direction  (for  example,  subtracting  39.7  gpm  from  30  gpm  gives  -9.7  ms'1  (reversed 
direction  of  the  gradient)  resulting  in  a  small  V  in  the  opposite  direction). 

2.3  Errors  Caused  By  Curved  Flow 

One  may  compute  the  difference  in  wind  speed  caused  by  the  flow  being  curved  instead 
of  straight  (along  a  great  circle).  For  this  comparison,  the  gradient  wind  for  moderate 

geostrophic  winds  (15  and  30  ms-1)  is  computed.  From  tables  in  List  (1984)  an  rf 
parameter  of  0.062  at  a  latitude  of  45®  is  obtained  for  a  radius  of  curvature  (r)  of  about 
600  nir.  («  1110  km).  For  Vg  =  15  ms'1, 

V  =  IU.6ffA~1ftm*~1  -1  +  (l  + - 4(29.1nmftj - .)*  =  26.1  An  «  13.5  ms~x  . 

1,  360<kr10.062nmsr,J 

Also,  for  Vg  =  30  ms1 , 

V  =  47.9  fat  «  24.7  ms1 

The  two  differences  (V  -  computed  above  are  1.5  and  5.3  ms~l . 

2.4  Errors  Caused  by  Mean  Temperature  Error 

Another  possible  cause  of  error  is  an  incorrect  value  of  T  but  no  difference  in  T  error 
between  the  two  soundings  (same  T  error  at  both  locations).  The  constant  level  surface 

version  of  equation  (1)  is  used;  that  is,  V  *  where  the  assumption  is  that 

fP 

Ap  «  20  mbar ,  tn  =  1000  km,  p  »  700  mbar,  tat  -  45® 

For  T  =  270  K, 

V  .  (287  K-Ws-*)  270  K(20  mbar\  m  2J  5  ^  > 

*  l.OSflO^'1  700  mbar\  10*  m  ) 
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and  for  T  =  267  K, 


Vg  *  21.3  ms1  . 

The  difference  in  these  moderate  values  of  Vg  is  about  0.2  ms'1 .  The  likely  error  of 
a  few  tenths  of  a  meter  per  second  for  moderate  V  is  small  relative  to  the  other 
possible  errors. 

2.5  Errors  From  Ageostrophic  Deviations 

A  significant  potential  source  of  error  arises  from  not  considering  certain  deviations  from 
the  geostrophic  wind,  which  may  be  roughly  approximated  by  the  horizontal  ageostrophic 
wind  components.  Equation  (4)  may  be  used  to  approximate  the  ageostrophic  component 
arising  from  the  pressure  tendency. 


Here  a  pressure  tendency  normal  to  the  wind  of  -1  hPa  over  1  h  near  45°  latitude  is 
considered.  One  hPa  =  10  jig  m'ls'2  and  the  air  density  (p)  =  1  kg  m'3 . 


ir  -  v  \  » - 

'  ‘  1  ig  m 


\ 


If  the  tangential  component  was  the  same  size,  the  rough  magnitude  of  the  total  pressure 
tendency  term  would  be  about  1.7  ms'1 . 


Equation  (5)  expresses  the  ageostrophic  component  arising  from  a  change  in  the  spacing 
of  contours  on  a  constant  pressure  surface  (or  isobars  on  a  constant  height  surface).  A 

change  of  geostrophic  wind  of  2  ms'1  over  a  distance  of  200  km  is  not  unreasonable. 
For  geostrophic  wind  speeds  of  15  ms-1  and  30  ms'1 ,  at  45*  latitude, 


(V  * 


15  ms'1 

ix^io-V* 


2  ms 


2(10 > 


*  1-5  ms’1 


and  (V  -  VJt  *  2.9  ms1 


Higher  values  of  Vg  or  sVg  can  occur,  leading  to  higher  values  of  this  "contour  spacing" 
component.  When  both  components  are  combined,  the  ageostrophic  part  of  the  "real" 
horizontal  wind  may  exceed  5  ms"1  (or  <  1  ms~l  if  the  two  terms  tend  to  cancel  or  Vg 
is  small). 
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2.6  Apparent  Error  From  Direction  Deviations 


So  far  the  error  caused  by  incorrect  orientation  of  the  gradient  has  not  been  explicitly 
considered.  Jedlovec  (1985)  reported  on  differences  between  gradients  derived  from  data 
of  the  visible  infrared  spin  scan  radiometer  (VISSR)  atmospheric  sounder  (VAS)  carried 
on  the  geostationary  operational  environmental  satellite  (GOES)  and  those  computed 
using  data  from  a  special  dense  net  of  rawinsondes  (50-km  spacing).  The  differences 
ranged  from  a  few  degrees  (almost  parallel)  to  about  90°  (nearly  perpendicular).  In  the 
absence  of  other  factors,  an  error  in  direction  of  the  gradient  (or  contours)  of  around  30° 
could  lead  to  a  difference  in  V  along  the  expected  wind  direction  of  about  2.0  and  4.0 

ms"1  for  V  values  of  15  and  30  ms'1,  respectively.  A  realistic  high  deviation  in 
gradient  near  45°  could  produce  Vg  differences  of  about  4.4  and  8.8  ms'1  for  the 
aforementioned  V  values.  While  even  good  quality  data  from  rawinsondes  may  contain 
wind  velocity  errors,  the  deviations  in  direction  reported  by  Jedlovec  (1985)  may  reflect 
real  differences  from  the  true  direction  of  the  gradient  (or  contours). 

2.7  Net  Error  Values 

Table  3  lists  some  of  the  potential  causes  of  error  in  satellite-derived  estimates  of  wind 
speed  that  use  the  geostrophic  wind,  along  with  possible  magnitudes  of  those  errors.  A 

moderate  wind  speed  is  assumed  (10  to  30  ms* ). 

TABLE  3.  POSSIBLE  MAGNITUDES  OF  ERROR  IN  WIND  SPEED  ARISING  FROM 
LISTED  CAUSES  FOR  SATELLITE-DERIVED  GEOSTROPHIC  WINDS. 


Cause 

Resultant  Wind  Speed  Error 
(ms'*) 

1 

Incorrect  T 

0.1  to  0.4 

2 

Gradient  instead  of  Geostrophic  Flow 

l  to  5 

3 

Computed  contour  gradient  incorrect 
(varying  error  in  T  sounding)* 

5  to  20 

.  .  -  --  .  .  . 

1 

Ageostrophic  component 

Pressure  tendency 

Contour  spacing  change 

0.5  to  2 

1  to  5 

5 

Incorrect  direction  of  gradient  (height 
or  pressure) 

Near  0  to  10 

(Component  in  expected  direction  of  wind) 

*For  the  T  sounding  error  (item  3)  a  layer  of  at  least  300  hPa  was  assumed  (for  example,  near 
sea  level  surface  to  about  3  km  AGL,  the  region  of  largest  satellite  temperature  error  over  land). 
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These  values  may  tend  to  augment  or  cancel  one  another.  For  example,  items  3  and  4 
may  result  in  errors  of  over  20  ms'1  or  only  1  or  2  ms'1 .  Miers  et  al.  (1992)  and  their 

references  report  satellite  geostrophic  wind  errors  from  ±4  to  ±14  ms'1 ,  generally 
compared  with  wind  speeds  measured  by  rawinsondes.  Rawinsonde  data  also  may  have 

errors,  ranging  from  <  1  ms;'1  for  high  quality  equipment  to  >  2.5  ms~l  for  older 
(now  mostly  obsolete)  equipment  (Fisher  et  al.  1987).  Reasonable  "middle  of  the  road" 
errors  (1  to  4  of  table  3)  from  this  analysis  may  lead  to  a  "total  error  of  about  15  to  20 

ms'1  where  the  individual  errors  sum  together  or  perhaps  only  2  or  3  ms'1  where  they 
largely  cancel. 

Deviations  in  direction  of  the  gradient  or  contours  from  the  actual  direction  (Miers  et  al. 

1992,  Jedlovec  1985)  could  result  in  apparent  errors  of  up  to  10  ms'1 .  This  "error" 
would  be  important  for  computing  the  wind  speed  for  a  specified  direction.  An  example 
would  be  the  cross  or  along  trajectory  wind  speed  for  an  artillery  application.  The 
apparent  error  (but  a  real  component  error)  could  either  increase  or  decrease  the  "true" 

error.  For  example,  an  error  in  wind  speed  of  10  ms'1  combined  with  a  deviation  in 
gradient  direction  could  result  in  a  down  range  error  of  approximately  15  or  5  ms'1 . 

The  net  errors  shown  in  table  3  represent  a  range  of  "typical"  values  for  moderate  wind 
speeds  (10  to  30  ms'1).  However,  the  reader  may  compute  other  values  of  potential 
errors  for  greater  or  lesser  wind  speeds,  as  well  as  for  other  parameters  (for  example, 

different  radii  of  curvature,  or  differences  in  T  error  less  than  0.5  K  or  greater  than 
2  K).  The  reader  may  elect  to  calculate  new  potential  errors  for  one  or  more  of  the 
causes  of  table  3  to  arrive  at  additional  net  error  values.  Appendix  A  con  dns  a 
selection  of  tables  from  List  (1984)  that  should  suffice,  along  with  the  equations 
presented  herein,  for  nearly  all  additional  calculated  errors. 

3.  THERMAL  WIND 

3.1  Basic  Derivations 

A  number  of  texts  provide  information  on  formulations  of  the  thermal  wind  (Haltiner  anc 
M  min  1957,  Holton  1979).  The  form  derived  in  Holton  (1979)  is  used  here.  Given  the 
horizontal  components  of  the  geostrophic  wind  in  geopotential  form, 


IQ 

fdy 


and  v. 


IQ 

fdx 


(7) 


and 


8<J>  1  RT 

—i  -T  «S  —  ...  SS  —  ■ 

dp  p  p 
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When  equation  (7)  is  differentiated  with  respect  to  p  and  equation  (8)  is  applied,  the 
results  are 


fry  _  fry  _  ... 

atop  /la*/, 

p*£  .  *».  .  ^ 

dbip  f\dy)p 
dV 

or  in  vector  form  — -  * - fcr(v7Y 

dlnp  f  > 

and  Xr  -  ^0>,)  -  K,0>„)  -  -  |  f'U  x  vT>dbp 

where  T  =  temperature,  p  =  pressure,  p  =  density,  R  =  gas  constant  for  dry  air,4> 
=  geopotential,  and  Vf  =  geostrophic  wind,  f  is  the  coriolis  parameter  (=  2  Q  sin8 

where  Q  is  the  angular  velocity  of  the  earth  at  the  poles  and  Bis  latitude).  In 
component  form 


where  uT  *  ujp j)  -  ujpj  and  vT  **  vjp j)  ■  are  the  thermal  wind  components 
defined  as  the  differences  in  the  geostrophic  components  for  the  upper  p  surface,  denoted 
with  subscript  1,  and  the  lower  p  surface,  denoted  by  subscript  o.  T  is  the  mean 
temperature  for  the  layer  from  p0  to  pl ,  and  x  and  y  are  distances  in  the  east-west  and 
north-south  directions,  respectively.  An  alternative  method  is  to  use  the  geopotential 
form,  or  the  equivalent  "height"  form. 

. .  i  q*.  -  ♦.) . .  -«.) 

r  ^ 

T  f  dx  f  dx 


Since  this  paper  is  concerned  with  errors  in  satellite  temperature  soundings,  equations 
(9)  and  (10)  will  be  used. 

3.2  Errors  From  Incorrect  Mean  Temperatures 

Equations  (9)  and  (10)  may  be  approximated  in  finite  difference  form  to  a  sufficient 
accuracy. 


Taking  differences  in  errors  in  T  between  soundings  200  km  apart  of  0.5,  1,  and  2  K, 
one  can  compute  potential  discrepancies  in  thermal  wind  for  given  layers  of  atmosphere. 
The  700  to  300  hPa  layer  will  be  considered  as  a  whole  and  subdivided  into  100  hPa 

sublayers.  Isotherms  orientation  will  be  assumed  to  be  NE-SW  (45°)  and  latitude  to  be 
43°.  For  the  entire  layer,  for  a  2-K  error  in  T 

um-  2*ZJT1|w2*~2f.,  1415  K  Vo.8473)  -  -  28.7(0.8473) 

HT4*1 
*  -  24.3  iw-1 
vr  «  24.3  ms'1 

The  magnitude  (wind  speed)  from  the  standard  formula  is 

VT  •  \i 4  +  *  34.4  ms1 

For  a  1-K  error  in  T , 
uT  -  12.15  mi'1 
vT  ■  12.15  ms~l 
VT  -  17.2  Mr* 
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For  0.5-K  error  in  T, 


VT  m  8.6  ms'1 

Table  4  shows  the  values  for  individual  sublayers  (I)  and  the  cumulative  error  (sum  of 
errors  from  lowest  to  given  layer,  C).  Slight  differences  between  last  cumulative  values 
and  values  for  entire  700  to  300  hPa  layer  appear  because  of  round  off  of  table  values. 

TABLE  4.  THERMAL  WIND  VALUES  (wj1)  FOR  LISTED  LAYERS  AND 

DIFFERENCES  IN  ERRORS  IN  MEAN  TEMPERATURE  (7) 
BETWEEN  SOUNDINGS  200  KM  APART. 


T 

Difference 

Layers  (hPa) 

700-600 

600-500 

500-400 

400-300 

I  C 

I  C 

I  C 

I  C 

0.5 

1.6 

1.6 

1.8 

2.3 

m 

2.9 

8.6 

1 

3.1 

3.1 

mi 

6.8 

11.3 

5.8 

17.1 

2 

6.3 

6.3 

MEM 

13.7 

9.1 

22.8 

11.7 

34.5 

3.3  Net  Potential  Error 

The  results  shown  in  table  4  do  not  include  possible  errors  in  the  baseline  wind  field  (for 
example,  wind  at  700  hPa),  errors  arising  from  not  considering  curvature,  and  ageo- 
strophic  components.  Use  of  measured  wind  as  the  baseline  value  will  reduce  the  total 
possible  error.  However,  since  the  thermal  wind  itself  makes  use  of  the  geostropic 
approximation,  the  other  sources  of  error  may  have  an  increasingly  greater  effect  as 
height  from  the  baseline  level  increases. 

As  a  first  estimate  of  potential  errors,  one  can  use  the  values  of  table  3,  with  the  range 
of  values  from  table  4  replacing  item  3  (varying  error  in  T  sounding  causing  incorrect 
contour  gradient).  Table  5  repeats  part  of  table  3  along  with  new  results  from  table  4 
for  a  700  to  300  hPa  layer.  The  values  of  item  1  of  table  5  would  be  smaller  for  thinner 

layers  (for  example,  a  range  of  <  1  ms'1  to  around  12  ms~l  for  differences  in  T  error 
of  a  few  tenths  of  a  degree  to  about  2  K  for  a  100  hPa  thick  layer). 

These  errors  may  tend  to  augment  or  cancel  one  another.  For  example,  items  1  and  3 

may  result  in  errors  of  over  35  ms-1 ,  or  only  1  or  2  ms* .  For  thinner  layers  (for 
example,  700  to  500  hPa)  the  value  of  item  1  would  be  much  smaller  (table  4),  reducing 
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the  overall  potential  error.  However,  the  numbers  in  table  2  assume  a  "perfect" 
measured  wind  velocity  at  700  hPa.  The  wind  speed  error  of  measured  wind  (rawinsonde 

or  radar  profiler)  is  likely  to  have  a  value  near  1  to  2  ms'1  (Miers  et  al.  1992,  Fisher 
et  al.  1987).  This  error  also  may  augment  or  diminish  the  "thermal  wind"  error.  Some 
rough  idea  of  potential  error  may  be  obtained  by  assuming  "middle  of  the  road"  errors 

from  table  5  (items  1  to  3).  The  potential  total  error  may  reach  values  over  25ms"1 

where  the  individual  errors  sum  together,  or  perhaps  less  than  10  ms~l  where  they 
largely  cancel. 

TABLE  5.  POTENTIAL  MAGNITUDES  OF  ERROR  IN  THERMAL  WIND 

GIVEN  POSSIBLE  DIFFERENCES  IN  T  ERROR  BETWEEN 
SOUNDINGS  200  KM  APART.* 


Cause 

Potential  Wind  Speed  Error 
(ms'1) 

1 

T  gradient  of  layer  (700-300  hPa)  incorrect 

5  to  35 

2 

Gradient  instead  of  Geostrophic  (curved  flow) 

1  to  5 

3 

Ageostrophic  component 

Pressure  tendency 

Contour  spacing  change 

0.5  to  2 

1  to  5 

4 

Incorrect  direction  of  temperature  gradient 

Ne;  r  0  to  10 

*The  assumption  is  that  the  difference  in  T  error  does  not  exceed  2  K  for  the  700  to  300  hPa 
layer.  Baseline  wind  measured  at  700  hPa. 

Deviations  in  the  temperature  or  height  gradient  from  the  actual  direction  (Miers  et  al. 
1992,  Jedlovec  1985)  could  result  in  apparent  errors  (but  a  real  component  error)  of  up 

to  10  ms~l  in  geostrophic  wind.  These  same  deviations  could  cause  errors  in  estimates 
of  thermal  wind  of  about  the  same  size  for  thicker  layers  (for  example,  700  to  300  hPa). 
Since  thinner  layers  (for  example,  700  to  600  hPa)  normally  would  have  smaller  thermal 
wind  values,  the  error  from  deviations  from  the  true  gradient  direction  would  be  smaller, 

perhaps  less  than  3  ms'1 .  If  a  measured  baseline  wind  was  not  available,  the  error  could 
be  considerably  larger.  For  example,  for  a  700  to  300  mbar  layer  use  of  the  geostrophic 

wind  as  the  baseline  wind  could  lead  to  errors  of  over  40  ms-1  if  the  net  thermal  wind 
error  and  the  error  in  the  baseline  geostrophic  wind  tended  to  sum  together. 
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4.  CONCLUSION 


Errors  in  satellite  estimates  of  the  real  wind  that  use  the  geostrophic  approximation  or 
thermal  wind  may  reach  significant  values.  Net  differences  between  the  actual  and 

geostrophic  wind  from  a  few  to  over  20  ms~l  are  possible,  with  reasonable  middle 

values  of  10  or  12  ms~l .  The  apparent  wind  speed  error  caused  by  differences  in 
orientation  of  the  geopotential  height  gradient  may  increase  (or  decrease)  the  component 
of  the  geostrophic  wind  in  a  specific  direction,  resulting  in  an  augmentation  or  reduction 
in  the  effect  of  the  other  "real"  errors.  This  effect  may  be  especially  important  for 
artillery  or  aviation. 

Net  differences  between  the  real  wind  and  that  derived  from  the  thermal  wind  may  range 
from  a  few  to  over  30  ms~l  for  thick  layers  (for  example,  700  to  300  hPa),  with 
reasonable  middle  values  around  15  to  18  ms~l .  For  a  thin  layer  of  approximately  100 
hPa  thickness,  the  errors  may  range  from  1  to  over  10  jws-1  with  reasonable  middle 

values  near  3  to  6  ms-1 .  The  apparent  wind  speed  error  caused  by  differences  in 
temperature  gradient  orientation  may  increase  (or  decrease)  the  component  of  the 
estimated  wind  in  a  specific  direction.  As  with  the  geostrophic  wind,  this  apparent  error 
may  augment  or  reduce  the  effect  of  the  other  "real"  errors. 

The  potential  errors  described  in  this  report  may  occur  when  satellite  data  are  used 
alone,  without  any  additional  data  to  "tie”  down  the  appropriate  satellite  values.  Some 
improvement  should  be  possible  through  the  use  of  new  sources  of  data  (Miers  et  al. 
1992)  such  as  radar  wind  profilers  and  accompanying  radio  acoustic  sounding  system 
(RASS),  ground-based  radiometers,  and  sensors  on  unmanned  aerial  vehicles,  combined 
with  processing  techniques  currently  under  development.  In  addition,  new  satellite 
sensors  will  eliminate  many  of  the  difficulties  associated  with  current  instruments.  For 
example,  new  active  sensors  planned  for  early  in  the  next  century  (for  example,  the  laser 
atmospheric  wind  sounder)  will  measure  wind  velocity  directly  for  cloud-free  lines  of 
sight,  and  new  passive  sounders  (for  example,  atmospheric  infrared  sounder)  will 
improve  temperature  profiles.  In  the  interim  the  combination  of  additional  sources  of 
data,  new  processing  techniques,  and  satellite  sensors  coming  on  line  by  the  end  of  the 
1990’s  (for  example,  special  sensor  microwave  imager/sounder  on  Defense 
Meteorological  Satellite  Program  satellites,  and  the  advanced  microwave  sounding  unit-A 
on  National  Oceanic  and  Atmospheric  Administration  satellites)  may  provide  the  best 
solution.  Descriptions  of  these  sensors  may  be  found  in  NASA  (1991),  Patel  (1992),  and 
Swadley  and  Chandler  (1992). 


18 


REFERENCES 


Fisher,  E.,  F.  Brousaides,  E.  Keppel,  F.  J.  Schmidlin,  H.  Herring,  and  D.  Tolzene,  1987, 
Meteorological  Data  Error  Estimates,  Document  353-87,  Meteorology  Group,  Range 
Commander’s  Council,  White  Sands  Missile  Range,  NM. 

Haltiner,  G.,  and  F.  Martin,  1957,  Dynamical  and  Physical  Met.  McGraw-Hill  Book  Company, 
NY. 

Holton,  J.R.,  1979,  An  Introduction  to  Dynamic  Meteorology.  2nd  Edition,  Academic  Press, 
New  York,  NY. 

Jedlovec,  G.,  1985,  "An  Evaluation  and  Comparison  of  Vertical  Profile  Data  from  the  VISSR 
Atmospheric  Sounder  (VAS),"  J  Adams  Ocean  Tech.  2:559-581. 

Jursa,  A.  S.  (Ed),  1985,  Handbook  of  Geophysics  and  the  Space  Environment.  Air  Force 
Geophysics  Laboratory,  Han  scorn  AFB,  MA. 

List,  R.  (Ed),  1984,  Smithsonian  Met  Tables.  Smithsonian  Institute  Press,  Washington,  DC. 


Miers,  B.,  J.  Cogan,  and  R.  Szymber,  1992,  Review  of  Selected  Remote  Sensor  Measurements 
of  Temperature,  Wind,  and  .MQistu£gfc..agdLCpmpariwnjQ  Rawiasonds  Measurements,  tr- 
0315,  U.S.  Army  Research  Laboratory,  White  Sands  Missile  Range,  NM. 

National  Aeronautics  and  Space  Administration  (NASA),  1991,  EOS  Reference  Handbook. 
Goddard  Space  Flight  Center,  Greenbelt,  MD. 

Patel,  P.  K.,  1992,  "Advanced  Microwave  Sounding  Unit-A,"  Preprints  of  the  Sixth  Conference 
on  Satellite  Meteorology  and  Oceanography.  American  Meteorological  Society,  Atlanta, 
GA,  462-465. 

Swadley,  S.  D.,  and  J.  Chandler,  1992,  "The  Defense  Meteorological  Satellite  Program’s 
Special  Sensor  Microwave  Imager/Sounder  (SSMIS):  Hardware  and  Retrieval  Algorithms, " 
Preprints  of  the  Sixth  Conference  on  Satellite  Meteorology  and  Oceanography.  American 
Meteorological  Society,  Atlanta,  GA,  457-461. 


19 


APPENDIX  A.  TABLES  FOR  CALCULATING  GEOSTROPHIC  AND  GRADIENT 
WINDS 


This  appendix  presents  a  series  of  tables  from  List  that  may  be  used  as  an  aid  in  calculating 
geostrophic  and  gradient  winds.  Accompanying  explanations  are  included.  For  further 
details  see  List  (1984), 1  Haltiner  and  Martin  (1957), 2  and  Holton  (1979).3 


‘G.  Haltiner  and  F.  Martin,  1957,  Dynamical  and  Physical  Met.  McGraw-Hill  Book  Company, 
NY. 

2J.  R.  Holton,  1979,  An  Introduction  to  Dynamic  Meteorology.  2nd  Edition,  Academic  Press, 
New  York,  NY. 

3R.  List  (Ed),  1984,  Smithsonian  Met  Tables.  Smithsonian  Institute  Press,  Washington,  DC. 
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0E0STR0RH1C  WINS,  CONSTANT  PRESSURE  SURFACE 
100  Biopotential  meter  eortaura 


The  salar  equation  for  the  fecscrcpiuc  wind  oa  a  constant  pressure 

r.-V* 


44 


‘fir 


whirs  *  is  the  feopotenixal  is  a  constant -pressure  surface,  a  is  distance  measured  in  the 
hcrfcsr.nl  direciaa,  /  is  the  coriolis  parameter,  and  V,  is  the  component  0:  the  jeostrophic 
wind  r.cmil  to  tha  direction  is  which  n  is  measured. 

Cs  a  constant  pressure  surface  with  contours  draws  for  intervals  at  ICO  geopoteatiai 
meters  (gpm.)  this  reduces  to 
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(To  convert  loots  to  other  measures  of  sossd  tea  Table  35.) 
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Tabus  »7  tcoevcuuoeai 

ai09TR0l»HIC  WINO,  CONSTANT  IMS48URX  SURPACX 
100  geapotentisl  metsr  contours 

when  An  is  the  contour  spigot  measured  in  dejrees  of  latitude  (Le,  one  unit  of  An  has 
the  lenzth  oi  one  detree  o i  latitude  at  the  place  (or  which  the  contour  spacer  >*  measured) . 
Tide  37  fives  values  o!  V,  in  Vcxts  u  a  function  of  in  and  latitude  with  auxiliary 
cslusns  fnriar  equivalents  of  A*  in  kilometers,  statute  Biles,  and  nautical  miles.  Ii  the 
latter  are  measured  by  a  sap  scale  true  at  some  other  latitude  the  value  should  be  cor* 
recced  to  the  latitude  at  which  the  measurements  are  taken  (tec  Table  165). 

Since  the  gtostrophie  wind  is  inversely  proportional  to  the  contour  spacing  and  directly 
proportional  to  the  contour  interval  (At  gpm.),  values  of  V,  tor  1/10  of  the  indicated 
spacing  may  be  found  by  multiplying  the  tabular  values  by  10,  etc,  and  tor  contour 
intervals  that  are  multiples  or  subsulttples  of  100  gpm.  by  multiplying  the  tabular  values 
by  Ad/  ICO  (e.f,  (or  2C0  gpo.  contours  multiply  by  2,  for  t0  cm.  ccntoun  multiply  by 
i.  etc). 
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Tamlz  40 
GRADIENT  WIND 

The  equation  for  the  jndier.t  wind  speed  V  is  cgi  units  is, 

■where 

rm  “radius  of  curvarur*'”  o:  the  trajectory  »R  tan  •  (see  Table  166), 
f  s Coriolis  parameter  2u  sin  9  (m  a angular  velocity  oi  rotation  oi  the  earth, 
Os  latitude), 

V,  zx.  geostropfcic  wind  speed  (sc*  Tables  37-39). 

Equation  (1}  can  be  rewritten  in  the  following  fans 


lor  any  of  the  following  consistent  combinations  of  units: 

a.  V  and  V,  in  miles  per  hour,  r  in  statute  miles, 

b.  “  “  “  **  knots,  r  in  nautical  aulas, 

c  "  *  “  *  kil  onset  ere  per  hour,  r  m  kQometara. 

V  is  a  function  of  the  parameter  *f  and  of  the  gaosuophk  wtndjspeed  V,  only.  In  applying 
equations  (1)  and  (2)  the  fallowing  sign  convention  is  necessary:  for  cyclonic  curvature 
rf  ^  Q.  for  tstiardocie  ccrvarurt  rr  <  u. 

Table  SO  A  gives  values  of  the  parameter  rf  as  a  faction  of  latitude,  e.  and  r.  Table 
40  B  gives  values  of  V  for  cyclonic  curvature  and  Table  40  C  for  eztisyclocic  curvature 
as  a  function  of  the  parameter  rf  and  the  goostrcphic  wind  speed  V+ 

To  find  the  gradient  wind  speed  at  •  given  point.' 

1.  Determine  the  latitude  0  and  the  value  r  of  the  trajectory.  (Table  166  indicates 

a  method  icr  finding  r  on  e  polar  stenographic  mag  projection;  for  other 
projections  an  estimate  most  he  made) 

2.  From  Table  40  A  find  the  parameter  rf.  (This  parameter  is  linear  in  r  so  that 

values  for  other  radii  than  those  given  may  be  readily  determined  from  the  table, 
^tg,  if  r* 2300,  add  the  values  for  rm 2000  and  r»300.) 

3.  D* tenant  the  gvostrophtc  wind  speed  V,  at  the  given  point  (see  Tables  37*39). 

(V,  and  r  must  be  is  one  of  the  cunsi«Bt  combinations  of  units  given  shove.) 

4.  Enter  Table  40  B  (evdenie  cue)  or  Tibia  40  C  (anticydonic  case)  with  the 

arguments  rf  and  V  The  corresponding  tabular  value  is  the  gradient  wind  V 
in  the  same  units  as  V 

(canton teg) 
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Tail!  40  (CONTINUED) 
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TABLE  40.  (CONTINUED) 
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TABLE  166 


RADIUS  OF  CURVATURE  ON  A  POLAR  STEREOGRAPHIC  PROJECTION 


In  computing  Trad::"*  w-lnd  speeds  (Table  40)  and  in  other  problems  it  is  necessary  to 
determine  a  factor  r  which  depends  on  curvature  of  the  trajectory.  This  factor  arises  in 
taking  account  of  the  horizontal  component  of  the  centrifugal  force  acting  on  a  particle. 
The  problem  is  twofcld:  (1)  to  determine  the  trajectory  of  the  particle  on  a  map,  ar.d 
(2)  to  determine  the  required  value  ot  r  if  the  trajectory  on  the  map  is  known.  The  first 
problem  is  of  such  rarer:  that  it  cannot  be  treated  adequately  here.  (Non.— In  many 
cases  an  approximation  is  made  from  the  curvature  of  the  isobars  or  streamlines.)  The 
second  problem  has  been  solved  for  the  case  of  a  polar  stereo  graphic  projection,  since  on 
this  projection  a  "small  circle”  on  the  earth  projects  as  a  dreie  os  the  map.  Table  40 
provides  a  means  for  computing  the  desired  r  lor  trajectories  on  a  polar  stereo  graphic 
projection. 

Let  R  be  the  radius  of  the  earth,  /  the  true  radius  of  the  "small  circle”  on  which  the 
particle  is  assumed  to  b:  traveling  a:  a  given  instant  and  •  its  ar.guiar  radios  (as  seen 
from  the  center  of  the  earth).  Then  r1  sin  &  Since  we  are  concerned  with  the  hori¬ 
zontal  component  of  the  centrifugal  force,  the  effective  horizontal  radios  of  the  curvature 
required  is  the  gradient  wind  equation  is  given  by  r  ssri  sec  a  *  R  tan  «.  If  an  are  on  a 
map  representing  die  instantaneous  trajectory  of  a  particle  «  air  is  determined,  this  arc 
may  be  regarded  as  a  portion  of  a  “small  circle.* 

To  determine  r  for  a  given  arc  of  a  trajectory  oo  the  map: 

1.  Complete  the  circle  by  extending  the  arc  (a  set  of  circular  templates  will  prove 

very  useful). 

2.  Find  the  meridian  which  passes  through  the  center  of  this  drcle. 

3.  Determine  the  latitudes  <*  and  *>  of  the  points  where  this  meridian  intersects  the 

circle  (extend  the  meridian  across  the  pole  if  necessary). 

4  A.  If  the  drcle  found  in  step  1  dots  net  cornua  the  pole,  find  the  difference  between 
and  <»j  and  enter  part  A  of  the  table  with  this  difference  as  the  argument. 
The  corresponding  tabular  valoe  is  the  required  radius  r  in  statute  miles,  from 
the  formula  r  =  R  an  i($»  —  $i). 

4B.  If  the  drcle  found  in  step  1  contains  the  pole,  find  the  sum  (e,  +  *»)  and  enter 
part  B  of  the  table  with  this  sum  as  the  argument.  The  corresponding  tabular 
value  is  the  required  radius  r  in  statute  miles,  from  the  formula  r  as  R  tan  [90* 

— i(*»  +  **)J. 

( contained) 
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TABLE  166 


RADIUS  OF  CURVATURE  ON  A  POLAR  STEREOORaPHIC  PROJECTION 
A.  Circle  sot  inciudicg  pole. 
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70 
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80 
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B. 

Circle  inclndiog  pole. 

4i  +  4e 
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6728 
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536” 
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70 
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5347 
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80 
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4633 

4532 

4473 
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4243 
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90 
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3889 
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3753 
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3440 
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3320 
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110 

2771 
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130 
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1721 
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1097 
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170 

346 
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69 
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